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ABSTRACT 



The generation of ambient noise by physical processes dependent on shearing of the 
boundary layer under sea ice is investigated. Special attention is paid to the identifica- 
tion of individual noise-generating mechanisms and the assessment of their relative im- 
portance. 

Recent studies of Arctic ambient noise are reviewed with specific reference to results 
showing particularly good or poor correlation between ambient noise levels and ice 
movement or relative current. Potential noise-generating mechanisms are described and 
categorized according to their small-scale driving forces and expected noise character- 
istics. More detailed quasi-objective investigations are then used to establish the relative 
importance of each mechanism as a contributor to the overall under-ice noise spectrum. 

Flow, Mechanical mechanisms, involving ice sheet fracture as a result of wind and 
current-induced bending moments, are found to be unlikely contributors. Conversely, 
processes in which ice fragments in current-driven motion under the ice interact to cause 
bumping and grinding noises, appear to be of probable importance. Turbulent pressure 
fluctuations in the boundary layer under sea ice are shown to be of significance at low 
frequencies on a local scale. The role of resonant cavities in the under-surface of the ice 
does not appear, however, to be an important one. 
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I. INTRODUCTION. 



Recent strategic thinking has placed great emphasis on anti-submarine warfare in 
the ice-covered waters of the Arctic Ocean (Young, 19S6). This emphasis is the chief 
driving force behind ever increasing efforts to fully comprehend and model underwater 
acoustics in this region. A significant proportion of this ellort has been directed towards 
the study of ambient noise under the ice and in the marginal ice zone. The importance 
of being able to predict ambient noise stems from its variability and its influence on 
passive sonar performance. 

The effectiveness of any form of passive sonar, be it shipborne, towed, air dropped 
or mounted on the sea bed, is described by the passive sonar equation. From Urick 
( 1 9 S 3 ) 

SL — TL = XL -DJ + DT 

where 

SL is the source level: a measure of the sound emitted by the target 

TL is the transmission loss: the amount of attenuation experienced along the acoustic 
path from the target to the receiver 

XL is the noise level: the level of background and self noise obscuring the signal at 
the receiver 

DI is the directivity index: a measure of the directivity of the sonar equipment 

DT is the detection threshold: the minimum signal level that can be detected by an 
operator using the sonar equipment. 

It is the ability to accurately measure or predict the magnitude of the terms in this 
equation which permits the effective determination of passive sonar performance. A firm 
knowledge of the parameters which govern the passive sonar equation therefore allows 
an operator to plan both offensive and defensive operations with an invaluable degree 
of foresight. An example might be the ability to forecast a time period, or geographic 
location, offering poor passive sonar conditions in which a submarine might safely con- 
duct noisy evolutions while avoiding detection by hostile forces. Conversely, a unit en- 
gaged in offensive operations against submarine targets under the ice would benefit from 
an ability to predict a good acoustic environment. 
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The magnitude of the source level term is dependent on the sound energy emitted 
by the object of the passive detection effort. It is therefore governed by engineering, 
naval architectural and tactical considerations not directly related to the environment. 
Transmission loss depends entirely on environmental factors, being a function largely 
of sea water and sea bed characteristics. Temperature and salinity gradients, both hori- 
zontal and vertical, are particularly important. Although this term is undoubtablv im- 
portant in fixing effective sonar ranges, variability over moderate distances and time 
periods is relatively small in Arctic waters. This means that good climatological, or better 
still, recently measured, profiles of temperature and salinity are usually adequate for so- 
nar range prediction work over a broad area and for considerable time periods. 

The noise level term is generally of the most crucial importance in the overall sonar 
performance prediction effort in Arctic waters. It depends to a large degree on the en- 
vironment and is known to be highly variable both spatially and temporally. Noise level 
incorporates the effects of self noise, created by the hydrophone itself or by the platform 
on which it is mounted, and ambient noise from the surrounding water mass. Ambient 
noise from the environment is a variable factor depending on many environmental forc- 
ing functions. Its accurate prediction is crucial to the successful assessment of passive 
sonar performance. 

The last two terms of the equation, directivity index and detection threshold, are 
functions of equipment design, operating technique and in the case of the latter, opera- 
tor effectiveness. 

Thus it can be seen that ambient noise prediction is an essential part of planning 
for Arctic operations. An additional utility can also be found in weapons, sensor, and 
underwater communications systems design, where a grasp of typical ambient noise lev- 
els is important. 

Considerable effort has been directed towards the modelling and prediction of am- 
bient noise under ice. The problem has been tackled in a variety of ways including em- 
pirical approaches, attempts to predict noise levels through the summation of emissions 
from sources and by consideration of energy budgets (Oard, 1987). Pritchard (1988) 
emphasised the need to investigate and understand the physics behind noise-generating 
processes on a local scale before modelling them and summing their effects at a partic- 
ular location. The overall predicted noise would then be the sum of the noise intensities 
from many locations, generated by a variety of mechanisms and subject to transmission 
losses depending on receiver depth, source depth, range to the source and frequency. The 
resulting Arctic ambient noise model would be driven by outputs from existing 



meteorological and ice dynamics models. A categorization scheme for noise-generating 
mechanisms was proposed, by Pritchard, based on forcing factors. Mechanisms impor- 
tant under compact ice conditions are separated into three areas 

1. ridging effects which are driven by large scale ice stress and deformation 

2. micro-cracking effects which are driven by local ice stress 

3. effects dependent on the shearing of the boundary layer of the ocean under the ice. 

It is the purpose of this study to identify and investigate any mechanisms which may 
fall into the last category. Such potential noise-generating mechanisms are catalogued 
as exhaustively as possible and categorized according to their small scale driving forces 
and expected noise characteristics. Investigations are then conducted with a view to 
establishing their importance, if any, in the overall generation of ambient noise under 
ice. 

Many of the physical processes involved are highly complex, and some are poorly 
understood or the subject of current research. For this reason the investigations are to 
some degree subjective. An attempt has been made to model the physics in as elementary 
a fashion as possible and generally to draw inferences from 'order of magnitude' results. 

Chapter II is a review of recent Arctic ambient noise research which considers the 
effects of boundary layer shear. The purpose of this brief review is twofold. Firstly to 
explore the fundamental concept that noise is caused by shearing of the boundary layer 
under ice. and secondly, to identify individual mechanisms which have been suggested 
as being responsible for noise. Evidence is sought to establish whether shearing in the 
oceanic boundary layer, in the form of ice movement through the water or water move- 
ment under stationary ice. is important as a forcing function for noise generation. 

Chapter 111 outlines potential noise generating mechanisms and a proposed cat- 
egorization scheme. The individual mechanisms are then investigated in Chapter IV with 
conclusions and a summary presented in Chapter V. 
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II. REVIEW OF CURRENT LITERATURE. 



Previously published studies of under-ice ambient noise are fairly numerous and in- 
vestigate a broad range of different forcing functions. Forcing by thermal effects were 
investigated by Ganton and Milne (1965) and Milne (1972) while forcing by pressure 
ridging and ice mechanics was studied by Pritchard (1984) and by Buck and Wilson 
(1986). Various other atmospheric forcing functions have also been studied by Payne 
(1964). Ganton and Milne (1965) and Greene and Buck (1978). 

Three recent papers, Makris and Dyer (1986), Lewis and Dcnner (1988) and Waddell 
and Farmer (1988). consider shearing of the oceanic boundary layer in the form of ice 
movement or flow of water under ice as a forcing function. These papers are therefore 
reviewed briefly with the objective of highlighting results which tend to confirm or deny 
the importance of boundary layer shearing under ice as a forcing function for noise 
generation. 

Makris and Dyer (19S6) analysed low frequency ambient noise measurements and 
comprehensive environmental data from the 1982 TRAM IV expedition to the central 
Arctic (S3°A’ 20°£). The octave band 10 - 20 Hz was taken as being representative of 
noise in the low frequency range and compared with temperature and various applied 
stress components. The underlying hypothesis was that:- 

Xoise is created by ice fracture mechanisms proportional to the state of stress in the 
ice as induced by environmental loading. 

Cross-correlation coefficients between the low frequency noise records and shear 
stresses due to the wind and current were found to be particularly high. 0.84. in both 
cases for 9.9-day records and 0.71 and 0.74. respectively, for 23.7-day records. The cor- 
responding coefficient for tensile stress due to cooling was 0.15. Composite measure- 
ments. namely ice stress on the ice sheet's vertical section and the stress moment acting 
about the ice sheet's central horizontal plane, were also highly correlated with ambient 
noise, correlation coefficients being 0.81 and 0.87, respectively, for 9.9 day-records. 

The general conclusion of the work was that low frequency pack ice noise correlates 
best with the moment due to opposing wind and current stresses acting on the ice and 
worst with air temperature. 
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Lewis and Denner (1988) studied the correlation between noise at 10, 32 and 1000 
Hz and a number of parameters describing temperature and ice kinematics. The data 
analysed were obtained during the AIDJEX project in the Beaufort Sea during 
1975-1976 and was exceptionally comprehensive. High quality positional data were 
available from some 40 satellite-tracked drifters which allowed the accurate determi- 
nation of ice translational speed and differential kinematic parameters. Correlation co- 
efficients were calculated for ambient noise levels at the above frequencies and the 
following parameters: 

L\ the translation rate of the ice 

L' 2 , proportional to the kinetic energy of the ice 

D, the divergence of the ice 

|D| 

V, the vorticity (rotation rate) of the ice 
|V| 

X, the normal deformation rate of the ice 

S. the shear deformation rate of the ice 

F, the temperature difference between the air and the ice 

and 

(.V 2 + S 2 ) 1 - 2 , the total deformation of the ice. 

The results of these calculations are reproduced as Figure 1. The correlation ob- 
tained throughout between ambient noise and U, the simple translation rate of the ice 
measured by satellite tracking, is quite striking. In virtually all cases the correlation was 
the maximum observed for any of the parameters above. U correlated particularly well 
with 32 Hz noise in the summer (August 1975). This was attributed to the ice "rushing 
through the water". Another observation was that, in general, U correlated better than 
U-, except at 1000 Hz. This suggests that pure ice motion is more important in produc- 
ing noise than the kinetic energy of the floes, except at higher frequencies where collision 
events are probably more significant. 
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Figure 1. Linear correlation coefficients between ambient noise and ice kinematic 
parameters (from Lewis and Denner, 1988). 



Overall this work seems to emphasize the importance of ice movement in the gen- 
eration of ambient noise. Although no particular mechanisms are mentioned, the role 
of shearing in the oceanic boundary layer appears to be of probable significance. 
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Waddell and Farmer (1988) describe the noise resulting from the break-up of land 
fast ice in Amundsen Gulf in the Canadian Archipelago. Although the study was in fact 
motivated by the possibility that ice conditions could be monitored acoustically, the 
conclusions are of some interest here. 

Data collected over a 127-day period from April to August 1986 included ambient 
noise over the frequency range 300 Hz to 14.5 kHz, current measurements and 
meteorological observations. Currents in the area are strongly influenced by the tides. 
This is shown in Figure 2 which also depicts time series of longwave infra-red, ambient 
noise and air temperature for a 10-day period in June. 
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Figure 2. Time series of ambient noise and various environmental parameters 
(Waddell and Farmer, 1988). 
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Initially, when the ice was consolidated and land fast, there was no obvious corre- 
lation between the current meter observations and ambient noise. A more marked cor- 
respondence was noted between noise and atmospheric factors. A strong northerly 
current, developing in July, however, coincided with a sharp increase in noise levels as 
the ice near the recording instruments Jpegan to break up (see Figure 3) L 




Figure 3. Time series showing the increase in ambient noise levels with increased 
current and ice breakup (Waddell and Farmer, 1988). 
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These results seem to indicate that when ice is landfast and static, a moderate cur- 
rent under it does not create discernable noise in the range 300Hz to 14.5 kHz. When 
the ice cover is in motion, however, significant correlations can be found between cur- 
rent speeds and high frequency noise. This is probably be due to collisions and cracking 
during ice breakup. 

As a general summary, the studies done by Makris and Dyer (19S6) support the idea 
that low frequency (10-20 Hz) noise results from opposing wind and current stresses 
acting on the ice. The results of Lewis and Denner (19SS) show strong correlation be- 
tween noise and ice translational speed which is proportional to the relative current un- 
der the ice. Waddell and Farmer showed that little noise results from relative current 
under stationary ice. 
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III. POTENTIAL NOISE GENERATING MECHANISMS. 



Ambient noise due to shearing of the boundary layer under sea ice could be caused 
by a number of different mechanisms. Little or nothing has been written about the de- 
tailed processes involved and in some cases the physical processes are sufficiently com- 
plex to prohibit meaningful analytical investigation. 

A number of potential noise-generating mechanisms are put forward here. The list 
is not intended to be exhaustive but most mechanisms which appear to have any possi- 
bility of being relevant to real world conditions are included. Some mechanisms appear 
intuitively more important than others while some may seem a little obscure. All are are 
included in this chapter and will be investigated to a greater or lesser degree in Chapter 
IV. 

A. A CATEGORIZATION SCHEME. 

Two broad categories of potential noise generating mechanisms have been identi- 
fied. Firstly. Flow/ Mechanical mechanisms which involve mechanical interaction between 
the sea and the ice cover and produce sound through cracking or impactive events. 
Secondly. Flow mechanisms which are purely a function of fluid flow in the presence of 
the solid ice boundary. 

'fhe potential mechanisms arc therefore divided and ordered in this study according 
to these criteria. Figure 4 summarizes the mechanisms that have been considered and 
their relative positions in this scheme of organisation. Figure 5 is a schematic summary 
showing the partition of Flow Mechanical and Flow mechanisms. 
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Figure 4. An overview of potential noise-producing mechanisms. 
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B. FLOW/MECHANICAL MECHANISMS. 

1. Cracking Noise Due to Wind/Current-Induced Moments. 

The concept of opposing wind and current stresses inducing moments in the 
Arctic ice cover was originated by Makris and Dyer (19S6). Although no detailed de- 
scription of the mechanism itself was presented, the high correlation obtained in the 
FRAM IV data, between noise and the stress moment acting about the ice sheet's hori- 
zontal plane, indicates the potential importance of this mechanism. 

In an attempt to investigate the detailed mechanics of the processes involved, 
the wind and current stresses are taken to have most effect on features with some ver- 
tical extent. Pressure ridges or ice hummocks present obvious points where winds and 
currents could act to set up bending moments in the ice sheet which might possibly lead 
to ice fracture and therefore noise. The magnitudes of such bending moments are a 
function not only of the relative velocities of the wind and under-ice current, but also 
of the size and shape of the ridge and of the thickness and general nature of the ice sheet 
itself. 

Maximum opposing wind and current velocities are most likely to occur when 
the ice is in a low-stress state with some degree of wind driven motion, Ice in a static 
state is less likely to be subjected to strong induced moments as the geostrophic and tidal 
currents in the open Arctic tend to be relatively weak. Wind and current stresses will 
also be greater in areas of highly deformed ice cover as drag there is more pronounced. 
Fracture as a result of induced bending moments may not ncccessarily be more likely, 
however, as significant deformation is often coincident with greater ice thickness, lower 
salinity and therefore greater strength. 

2. Ice Fragment Bumping. 

The importance of ice collisions in the production of ambient noise has been 
recognized in several works (Diachok and Winokur ,1974: Shepard, 1979) and is an es- 
tablished mechanism when considered on a large scale, i.e., floe-floe interactions. The 
presence and movement of significant amounts of ice debris under sheet ice is also doc- 
umented (Zubov, 1943; Buck and Wilson, 1986). A potential noise generation process 
therefore exists in any pack ice covered region where deformation events create quanti- 
ties of loose fragmented material which can be moved around by a relative current. 
Noise from this source is likely to be aurally detectable and characterized by bumping 
or grinding sounds. 
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C. FLOW MECHANISMS. 



1. Noise From Turbulence. 

Turbulence is a recognized source of low frequency noise in the open ocean 
{Wenz. 1962: 1972). Although the processes by which it is radiated are exremely ineffi- 
cient and its far field effects are minimal, the direct effect on hydrophones within tur- 
bulent flows can be significant. This is because the pressure fluctuations impinge directly 
on the active surfaces of hydrophones which react by generating an electrical signal. 
Thus turbulent pressure fluctuations can cause a type of self-noise whose effects may 
be significant at very low frequencies. 

In the Arctic the situation is complicated by the presence of ice cover. The ice 
acts as a boundary of variable roughness past which water flows. In addition, the ice 
sheet itself may play a role in the radiation of turbulence noise to ranges beyond the 
actual source region. Turbulent self-noise has been recognised as a possible contaminant 
of Arctic ambient noise data by Makris and Dyer (19S6) and is included as a noise- 
generating mechanism in its own right. 

Flow-induced noise and flow-excited noise arc described briefly in Chapter IV 
as potential mechanisms by which the effects of turbulent pressure fluctuations are ra- 
diated to the far field. 

2. Resonant Cavities. 

The under-surface of floating ice in Arctic regions is. by its very nature, rough, 
deformed and pitted with cavities caused by melting, fracture or brine rejection. Such 
holes or cavities present a potential noise source when water is in motion across their 
openings. Flow may be the result of ice motion or currents under static ice and could 
excite resonant responses in these cavities. 

No known reference has been made to the possibility that resonant cavities 
might be responsible for under-ice noise but similar above-water features are certainly 
capable of producing noise in strong winds. Noise from hull cavity resonance is also well 
known in naval architecture {Blake, 19S6). 

The parameters affecting the frequency and intensity of sound from resonant 
cavities would be cavity size, shape and water flow speed. 
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IV. AN INVESTIGATION OF POTENTIAL NOISE MECHANISMS. 



A. FLOW/MECHANICAL MECHANISMS. 

1. Wind/Current-Induced Moments at Pressure Ridges. 

The action of a current relative to an ice Field on submerged pressure ridge keels 
and the action of the wind on their above water sails gives rise to a potential 
flow mechanical noise-generation mechanism. The moment set up by these forces may 
be enough to cause cracking and cracking noise in the ice sheet close to the ridge fea- 
tures. 

The moment M k due to water pressure on the ridge keels is likely to be rein- 
forced by wind action on corresponding sail formations causing a moment M s {see 
Figure 6). This is particularly so when internal stresses within the local ice field are low 
and motion is largely due to wind-induced free drift. The applied moment M c is therefore 
a function of the wind velocity V w . the relative current velocity V e and the pressure ridge 
geometry. 3/„ is the sum of the moment caused by the relative current on the keel 
and M, the moment on the sail due to the wind. 

This functional relationship is investigated below by assuming an idealized 
pressure ridge model after Wright et al. (1978) and Paquette and Bourke (1988). The 
geometry of this model is illustrated in Figure 7. 

It is assumed that 

1. the system is in local isostatic equilibrium when V w = F c = 0 

2. the ridge is symmetrical about point O as shown in Figure 7, with dimensions 
governed by local isostatic balance 

3. the density of the ice is homogeneous throughout with an equal percentage of void 
in the sail and keel structures 

4. the surface of the ice sheet, sail and keel are perfectly smooth 

5. the wind velocity v H (z) and relative current velocity v c (z) vary linearly with height 
and depth as shown in Figure 8. 
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a. Ridge Model Geometry. 

From the assumption of isostatic equilibrium 



Z 2 h s lV sPi = (h k lV s + V2h k (lV k -lV s ))(p w -Pi) 



where 

p w is the water density 
Pi is the ice density 



( 1 ) 
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From Figure 7 



If 






Jh_ 

tan y. k 



and 



»> 



tan a f 



so 



2/?i 



JF, = t 

tan 



+ »'c 



Substituting (2) and (3) into (1) and rearranging 

/ Pi \ _ , /?, tan a, 

/;: V ~ Pi ) ’ ^ tan 



( 2 ) 



(3) 



(4) 



Solving this quadratic in 




p k 



1+1 + 



/ p[ \ tan«j 
\ Pvt' - Pi ) tan a* 

( ^ 

\ Pw Pi J 



1/2 



(5) 



Now taking typical values for />, and p„ as 920 kg nr 3 and 1026 kg in' 3 , respectively, and 
o. k = 30°, otj = 20° (from Wright et ah, 197S) 

h s 

-y- ^ 0.41 ( 6 ) 

"k 
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also from Fieure 7 



JK 



V - l k + hk + ~ : r tan a k 



or 



tan a i 



hi.' — hu + h s — — + 

Since the ice sheet itself can also be assumed isostatically balanced 
hPi = l k(Pw - Pi) 
and since 

‘ = h+ (ft 



( 7 ) 



( 8 ) 



't-'jt* 0 - 9 ' 



( 9 ) 



and 



(Pw ~ Pi) 



h i p — 0 . 1 / 

rw 



( 10 ) 



Substituting (9) into (7) 



tan y. k p< 

h k , ^h k + h s -rrrz L + t-‘ 



Further substituting numerical values for a,_ a* and h k from (6) 

h k~ Ah s + l k 



( 11 ) 



or 

V- l-7/t* + ^ 



( 12 ) 
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From the geometry of the sail 
h s ' = h s + l 5 



and from (10) 



V = h s + 



iPw - Pi) 

Pw 



( 13 ) 



so 

h s 's*h s + 0.U (14) 

Now that functional relationships have been established between the vari- 



ous ridge dimensions, an equation for M„ the moment caused by wind action on the sail 
about point O can be found (see Figure 9). 
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b. The Moment Caused by W ind Action on the Sail. 

Using concepts from Zubov (1943) and the theory for a fluid jet acting on 
a flat plate (Hannah & Hillier. 1970) the wind velocity v* at height z above sea level can 
be resolved into two components: v, the component parallel to the surface of the ridge 
sail and v„ normal to this surface. 

V P = V W COS Oy (15) 

b, = sin cf. s (16) 

At the windward surface of the sail v„ causes a dynamic force. v p has no effect if the sur- 
face of the sail is perfectly smooth. The resulting force P is equal to the loss of momen- 
tum to the sail per second. So 

P = im w sin oq kg m s -2 



where 

m is the mass of fluid (air) incident on the sail at height z. per second. 

Then 

P = p 0 v 2 sin a s kg m s -2 ( 1 7) 

where 

p„ is the air density. 

The moment m s {:) about point O is given by the product of the length of the 
lever arm { and the component P e , of P. normal to { (see Figure 10). 



?■> 




Figure 10. Moment about O, the point of symmetry, due to the wind. 



From Figure 10 
P e = P n cos /? 



and 

/? = nj2 — 6 + a s 



so 



Pf= P sin((? — a^) 



and from (17) 

2 

P { = p a v w sin a s sin(0 — aj (1 8) 
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Now sin(0 — a,) can be found in terms of h,', €, a, and z by considering triangle OAB, (see 
Figure 11). 



B 




Figure 11. Triangle OAB. 



From Figure 11 

sin(# — a 5 ) = sin 6 cos a s — cos 6 sin a s 
or 

x cos a t — z sin a, 
sin(0 - a s ) = 

where 

V - ^ 

x = — 

tan a s 

so 



sin(0 — a s ) = 



h s ' cos 2 a s — z 
£ sin a s 



(19) 
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Substituting (19) into (18) 



and 



Pa' u(V cos\ ~ -) 
£ 



(20) 



m s( : ) = f P e = Pa v l( h s cos - z ) 

If the wind shear is linear, as outlined in the assumptions above 

,, 

v w (z) = 7 



so substituting (22) into (21) 



V 2 

/ \ ft Q w s \2 /it 2 \ 

m s (z) = — — (2 - t s ) (h s cos y. s - z) 



( 21 ) 



( 22 ) 



(23) 



The total moment M s on the sail about point can now be found by inte- 
grating with respect to z between t s and h‘ 







V 

(z - tf(h s ' cos 2 y. s - z)dz 



(24) 



or 






Pa* 



12 /:; 



(/:/ 4 (4cos\ 



3 ) + V \(8 



12 cos^) 4- h s ' 2 t 2 (6 — 12 cos'aj — 4 h s 'i* cos“ 



+ h) 
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From (13) 

V = !l s + h 
so 

J rl 

M s = (h;i 4 cosx -3) + ^/,(4 cos 2 a, - 4)) 

or 

M s = p a vW s C s (25) 

where C s is a form coefficient dependent on the geometry of the ridge. 

Q = -jT ( 4 cos2 °h “ 3 + (4 eos 2 a f - 4) j (26) 

I f 20° and 0. 1 1 

Q.(44-4-i-)x 10 - 3 (27) 

Figure 12 shows C s plotted against h s for various values of t. Note that C s ^ 0.04 for sail 
heights greater than one metre and is virtually independent of t. If the slope of the 
sail, is allowed to van - , an interesting effect is noticeable (see Figure 13). Where 
a, > 28° , C s reverses sign, i.e., the result is a moment in the opposite direction. This 
could be relevant when considering young, relatively unweathered ridge features where 
slopes may be greater than the 20° assumed above. 
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c. The Moment Caused by Current Action on the Keel. 

The moment M k caused by the relative current on the keel can be found in 



a similar way to that found for the sail. The relevant geometry is outlined in 
Figure 14. 
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By comparison with (21) 



m k = PP i = PMc( h k COS °-k - z ) 



Linear current shear is assumed, so 



l 'c( 2 ) — 



h 

r/2 



m k( z ) = U - l k )\ h k cos \ ~ 2 ) 



and 



M k = 



P_jj_ 

4< 



Ph 



0 - l k) 2 U l k COS 2 a k - z)dz 



(K + '*) 



This eives 



-'4 — j i 4/^4/ cos a k — 34c — 



From (12) h.'^Uh. +t k . so 



Fvjf , , 2, 



-'4 = — — ( W-S cos~a k - 3) + V„(4 cos‘a fc -4)) 



If a, ^ 30° and /, ^ 0.9/ 



.!/,< = ^y- 1 2/1,; -0.9V) 



or 






where C is a form coefficient for the keel given bv 



C b * [ 163 - 75 j- ) x 10 -3 



(28) 



(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 
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Figure 15 shows C k plotted against h k for various values of t. There is more dependence 
on i than in the case of C s and it should be noted that Q^O.15 for h k > \ m. Also 
comparing Figures 12 and 15. C k is generally about four times bigger than C, . When 
a k is increased, there is a point at about 44° where M k reverses sign as in the case of 3/ r 
at greater slopes (see Figure 16). Again this may be relevant when considering relatively 
new, unweathered, ridges. 
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The total moment M c on the ridge structure about point O is given by the 
sum of the moments caused by wind on the sail and relative current on the keel. 

M a = M s + M k (37) 

Figure 17 shows M„ with its constituent components M, and M h , as given by (25) and 
(35), plotted against h, with t = 1.5 m. Typical values for p a , V w p w and V e have been used, 
they are summarized below: - 
P„ = 1.3 kg nr 3 
p w = 1026 kg nr 3 
F„ = 10 m s* 1 
V c = 0.2 m s-' 

Tliese wind and current velocity values were chosen to represent the upper limits of ob- 
served measurements from the literature , i.e., a worst case has been chosen in order that 
the resulting moments will be the largest that could normally be expected to occur in 
nature. 
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t (metres) 



Figure 18 is a plot of M a against t and h , with V W =10 ms -1 and V c = 0.2 ms-' . The 
weak dependence on ice sheet thickness is apparent as is the order of magnitude of M c 
to be expected for these conditions , i.e., M c ~ 1 kN m per m ridge length. 




hs (metres) 

Vw = IQm/s, Vc = 0.2m/s 

Figure 18 . The total moment applied by wind and current at a pressure ridge as a 
function of sail height and ice sheet thickness (kNm per m ridge length). 
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The following points arise from the preceding calculation of the moment 
due to wind and relative current action on a pressure ridge. 

1. M k is the dominant component of M a . 

2. t, the thickness of the ice sheet has little effect on the absolute value of M c . 

3. For angles of slope a r and a* taken from Wright et al. (1978) M a is anticlockwise 
when the wind blows from the left and the current acts from the right (as shown in 
Figure 6). If steeper slopes are assumed, \t a could be reversed. Regardless of the 
sense of the resulting moment, however, the absolute value of M a at the ridge will 
be of much the same order. 

There are two considerations which have not been taken into account when 
calculating M s and M k which might lead to increased moments at the ridge. 

1. The surfaces of both the sail and keel have been assumed smooth. If a more real- 
istic approach were to be taken and some roughness allowed, drag could be ex- 
pected parallel to the sail and keel surfaces. 

2. Reduced pressure in the flow downstream of both the sail and keel would also have 
some effect on M s and M k . This would contribute in a clockwise sense for a wind 
blowing from the left and current from the right. 

Despite these possible inaccuracies and any others which might arise from 
the assumptions made, equations (25) to (28) and (35) to (37) allow a reasonable esti- 
mate of M a to be made for sensible values of ridge height and ice sheet thickness. 

2. The Effect Of An Applied Bending Moment On The Ice Sheet. 

In order for the moment .)/„. described above, to cause ambient noise under the 
ice. cracking must occur in the ice sheet. The horizontal distribution of bending moment 
.)/(. x) on the ice sheet must therefore be investigated in order to estimate the maximum 
moment that could be expected under normal circumstances. Two extreme cases are 
examined to determine whether the maximum expected bending moment on the ice sheet 
is likely to be sufficiently large to cause cracking. The first case examines the ridge as 
a relatively new feature when it is in an unconsolidated state. That is, it is assumed that 
the ridge structure consists of loose rubble and ice blocks and therefore contributes no 
strength to the system. The second case, representing an older pressure ridge, occurs 
when the ice debris comprising the sail and keel has frozen into a consolidated mass 
which together with the ice sheet near the ridge can be considered rigid. Figure 19 
compares, schematically, the two conditions outlined above. 
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Loose Material 



Unconsolidated Ridge 




Consolidated Ridge 



Figure 19. Consolidated and unconsolidated pressure ridges. 
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a. The Unconsolidated Pressure Ridge. 

In this case the ridge ice sheet system is assumed to behave in the same 
way as two adjoining, loaded, semi-infinite beams on an elastic foundation. The uncon- 
solidated mass of rubble and ice blocks forming the sail and keel interact with the wind 
and current as described previously to cause an applied moment M a and also a distrib- 
uted loading on the ice sheet. They do not however contribute in any way to the strength 
of the sheet. In order to calculate the expected bending moments on the ice sheet the 
section to the right of the ridge axis is treated as a semi-infinite beam with static loadings 
contributed by the sail and keel as shown in Figure 20. The total loading can be thought 
of as being due to three triangularly distributed loadings. That is, a downward load due 
to the weight of the ice in the sail (triangle ABC in Figure 20) and an upward loading 
due to the buoyancy of the keel (triangle DEM minus triangle FGH). So in the isostatic 



condition 

jy (. v ) = y ABC — }'deh + Xfgh (3 S ) 

where y is the downward displacement of the ice sheet, or 

=y<Jo^) -y { u P ) ( 39 ) 

Similarly 

= M a bc ~ ^h)EH + fgh H (l ) 

or 



M is the bending moment on the ice sheet (positive in the clockwise direction). 
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Equations for the displacement and bending moment on a semi-infinite 
beam subject to triangularly distributed loading have been developed from Hetenyi 
(1947) in the appendix. They are:- 

y(.v) = q 0 F{k, F,x) 0<x<F (42) 

where F( F,x) = (C iU . x , + 4J.(F — x) + B lx A l( — D lx (A xe +2 B u )) 

y(.v) = q 0 G{/.. f,x) x>F (43) 

where G(/., .v) = (C^i <f _ J ,| + B lt A )e - D lx {A )e +2 B i( )) and 

M{x) = - q 0 H{/, t?,x) 0<x (44) 

where //(/, F.x) = (A xU _ xl — D )x A )e — B )X {A )( + 2B, { )) 

{ is the base length of the triangular loading 

and <7 0 is the loading at the free end of the beam, i.e., at .v = 0 

see the appendix for definitions of 1// the characteristic length, k the foundation 
modulus and the functions A x . B x . C, and D x . 

The net downward loading q 0(d<!un) is found by considering triangles ABC and 
TGI I (see Figure 20). 

%(down) = h sPiS + (V - h k ~ l k)(Pw ~ Pi)g ( 45 ) 



and using (11) and (6) 

%(down) = h sS(Pi + 1 - ( <P* ~ Pi)) 

Taking p, ^ 920 kg nr 3 and p n ~ 1026 kg nr 3 
%(down) = 10670/?, kg m' 1 s' 2 

From Figure 20 

l V s h s 

^ (down) 2 _ tan sq “ 2,S/! * 



(46) 



(47) 



3S 



For the upward loading q 0{up) (due to triangle DEH) 

% (up) = (V - l k)(P w ~ Pi)S ( 4S ) 

% (up) = 4h s(Pw ~ Pi)S 

or, using the constant values given above 

% (u P )- 4l55h s kg m -1 s -2 (49) 



From Ficure 20 



IV k 

s — — — 

f (up) o 



h k 



tan y.fc 



+ 



IV, 

— - ~7/i 

-) ' n S 



(50) 



Substituting for q 0 and f in (42). (43) and (44), also taking A' = (/?„ — P,)^ — 1039 
kc nr 2 s~ 2 



y«tomp) = °- 942 r(/ - f (down)- x ) 


0 < x < € {dmtn) 


(51) 


y(Jown )( x ) = °- 942 G(/. f (downy x ) 


x — G (down) 


(52) 


2 ^(dow/j)( x ) k'5 //(/ , f (down)' •'") 


0 < X 


(53) 


3' ( ^)W = 0-1 44 F{/., f (up) ,x) 


o<x<y \ up) 


(54) 


y { u P )( x ) = 0.144 go.. y (upY x) 


■ X ^f(up) 


(55) 


•'WW = ~ 75 f(upy x ) 


0 < x 


(56) 
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Both <f (v) and F idonn) are related to h s by (47) and (50) and 7: is a function of i, so from (39) 



V/M = 0.942 F (down) (t. h s , x ) - 0.144 F (up) (i. h s , x) 
JV(-v) = 0.942 G {down) {t. h 5 , .v) - 0. 144 F (up) {t, h s , x) 
yfa) = 0- 942 G {dow4 1 ' h s- x ) - 0. 144 G (up) {i, h s , x) 
and from (41 ) 

Mfa) = — } S5 H (down) {i. h s ,x) + 75 H {up) {i, h s ,x) 



0 — x < ^ (down) (■ 57 ) 
^ (down) — x < ^(up) ( 58 ) 

x >F {up) (59) 



0 < .r (60) 



Figures 21 and 22 show the displacements of ice sheets of various thick- 
nesses when subjected to the loading of unconsolidated ridges with sail heights of 1.5m 
and 3m. It can be seen that the concentrated mass of the sail causes downward dis- 
placement at the ridge axis. The more distributed buoyancy force due to the keel causes 
upward displacement away from the ridge. This results in bending moment curves as 
shown by Figures 23 and 24. The value of the constant 7 was calculated using a typical 
value for E, the Young's Modulus for ice. E— 3 x 10° Pa. from Mellor (1986). 

The absolute value of the bending moment due to the isostatic loading on 
the ice sheet reaches a maximum value M mn some 10 to 20 m from the ridge axis. It is 
at this point, with or without the additional applied moment .\/ c due to the wind and 
relative current, that cracking is most likely to occur. Figure 25 shows .W m „ plotted 
against h s and i. 

Having established estimates for M a the applied moment and M m „ the 
maximum bending moment on the ice sheet due to isostatic loading, the two can be 
compared over a range of ridge heights and ice sheet thicknesses. From Figure IS and 
Figure 25 it is apparent that M c is of the order 0 to 1 k\ m per m ridge length for real- 
istic maximum wind and current speeds. M max , however, is of the order 20 to 2000 kN 
m per m ridge length, i.e., ,\/ max is roughly two orders of magnitude greater than M c for 
a given ridge size. 

Thus, it is concluded that if cracking occurs in the unconsolidated case, it 
will be due almost entirely to the isostatic loading on the ice sheet. The additional ap- 
plied moment M e , while not being sufficient to cause cracking by itself, could possibly 
act as a trigger mechanism for cracking where a cracking event due to isostatic loading 
was likely to occur anyway. 
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t = 1.5m 
t = 2.0ra 



Figure 21. Ice sheet displacement due to isostatic loading of an unconsolidated 



pressure ridge, height 1.5m. 
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In order to establish the likelihood of cracking due to isostatic loading 
alone, M maK is compared to the critical bending moment M cr „ at which cracking will oc- 
cur. M„„ is dependent on o, the tensile strength of the ice sheet, and its thickness t. A 
nominal value for the tensile strength of the ice sheet was taken, o, = 0.4 x 10 6 Pa from 
Mellor (1986). 

of 

Me* = ( 61 ) 

Figure 26 shows the ratio plotted against h, and t. From this graph it can be 

seen that 1 where t ~ h, . When M m JM er(r > 1 cracking due to isostatic load- 

ing is likely to occur. It should be noted that as ft, becomes much greater than t the result 
is a greater likelihood of cracking. 
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0 . 0-1 . * . » » 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

hs (metres) 

Figure 26. The ratio between maximum bending moment and critical bending mo- 
ment for an unconsolidated ridge under isostatic loading as a function of 
sail height and ice sheet thickness. 



44 



b. The Consolidated Pressure Ridge. 

In this case a completely rigid system is envisaged. It is assumed that the 
ridge structure has undergone repeated freezing events and is no longer active in that 
ridge building processes have ceased. The material comprising the sail and keel has fro- 
zen into a solid mass and contributes fully to the strength of the ice sheet. The solid ridge 
section is therefore treated as a rigid beam pivoting about point O (see Figure 27) while 
the adjoining ice sheets either side of the ridge are assumed to act as semi-infinite beams 
on elastic foundations. The boundary at section A-A, between the rigid portion (the sail, 
keel and an unknown length of ice sheet) and the rest of the ice sheet (which is consid- 
ered flexible), lies a distance L from the ridge axis. Equations for the displacement and 
slope of a semi-infinite beam on an elastic foundation are taken from Hetenvi (1947) and 
can be found in the appendix. 

Tor simplicity, symmetry about the ridge axis is assumed such that 

1. The distances L l and L R equal. 

2. y L andy*, the displacements of the ice sheet at distances L either side of the ridge 
axis, are equal and opposite. 

3. M l and X1 R . the bending moments on the ice sheet at distances L either side of the 
ridge axis, are equal and opposite. 

4. 0 L and Q R . the shear forces at distances L either side of the ridge axis, are equal and 
opposite. 

Thus 

L l = —L r = — L 

yL = -XR=y 
m l = -m r = m 

and 

Ql = ~Qr = Q 
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Figure 27. Geometry and general structure of the consolidated pressure ridge model. 
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Now summing the moments about point O 
S.U 0 = M b + \l L - \l R - LO l - LQ r - M a = 0 



or given the assumption of symmetry' above 



M b + 2M + 2LQ = 0 



( 62 ) 



M b is the moment caused by the immersion of ice to the left of O and emmersion of ice 
to the right. 



From the appendix, the displacement of a semi-infinite beam subjected to a shear force 
O and moment M at its free end is given by 



See the appendix for definitions of 1// the characteristic length, k the foundation 
modulus and the functions A x< B t C, and D x . Since the section A -A, distance L from the 
ridge axis, coincides with the begining of the flexible ice sheet and therefore the free end 
of the semi-infinite beam model, .v = 0 and D )x = C )x = 1. Therefore 



M b = x 2 k tan 6 dx = 2/3/<L 3 tan 6 



or 



M b ^ 2/3/tL 3 0 for small 6 



(63) 




(64) 
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The slope 6 of the semi-infinite beam section is also given in terms of Q and \f in the 
appendix. 



2 O) 2 



r * 3 



^ J A /.x + D;.x 



Again at the free end .y = 0 so A x , = D x , = 1 



20 /“ J \/> 3 '» / i 2 



Also, from Figure 27 



(65) 



tan 0 = y- or y = 0L for small 6 
Now subsituting (63) and (65) into (62) 
y L 3 / 2 (2.U/. - O) +2 m +2 LQ - M a = 0 

and from (64). (65) and (66) 

-U/( 1 + 2/L) 

® ( 1 + /. L) 



(66) 



(67) 



(6S) 



so (67) becomes 

M( y (f.L) 3 + 4(;.L) 2 + 4(/L) +2) -34(1 + (/L)) = 0 



The efiect of the applied moment on M, the bending moment at the interface region 
between the rigid pressure ridge section and the flexible ice sheet, is therefore given by 



_14 = U+(/T)) 

*' /fl ( y (/• T) 3 + 4(/. L) 2 + 4(T L) + 2) 



(69) 



Figures 2S and 29 show the ratio plotted for sail heights of 2 nr and 4 m with 

various ice sheet thicknesses. /. has been calculated using a Young's Modulus of 
E= 3 x 10 9 Pa and wind and current speeds of 10 ms' 1 and 0.2 ms' 1 have again been 
assumed to represent a worst case. 
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Figure 29. The ratio of applied moment to the maximum bending moment on the ice 
sheet for a consolidated pressure ridge, height 4m. 
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Note that (69) is not valid for L since the keel and/or the sail struc- 

tures significantly increase the rigidity there. It is assumed that cracking is unlikely 
within the ridge structure itself due to its increased thickness and strength. From 
Figure 29 the maximum bending moment on the ice sheet itself appears to occur where 

L = — 

2 ' 

The maximum bending moment M m „ is plotted for various values of sail 
height and ice sheet thickness in Figure 30. Note that M mtx is of the order 0 to 0.5 kN 
m per metre of ridge length. Comparing this with M cr „ , the bending moment required for 
cracking to occur as illustrated by Figure 31, the bending moment available is about 
three orders of magnitude smaller than that required for cracking to occur in the ice 
sheet. 



Thus, it is concluded that cracking noise is unlikely to occur due to 
wind/current induced moments applied to a consolidated pressure ridge. 




induced moments at a consolidated pressure ridge (kN m per in ridge 
length). 
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function of thickness (kN in per m ridge length). 



3. Conclusions : Cracking Noise Due to Wind/Current-Induced Moments. 

From the results of a) and b) above it is apparent that the generation of noise 
due to these processes is unlikely to be of any major significance, if it occurs at all in 
nature. In the case of the young, unconsolidated, pressure ridge it appears that the 
bending moments on the ice sheet due to the isostatic loading of the sail and keel 
structures far outweigh any additional bending moments caused by the effects of wind 
and current on the ridge. Any cracking which occurs is likely to be due to these loadings. 
An additional contributor to local ice sheet bending might also be large scale ice stress. 
This has not been considered here as it w'ould probably be small if the ice were in wind 
driven motion. There remains the remote possibility that the cracking of heavily loaded 
ice sheets could be triggered by the onset of wind/current-induced effects. It seems un- 
likely however that this could account for any significant correlation between wind, 
current and ambient noise in field measurements. The results of the investigation into 
the more mature consolidated pressure ridge reveals even less probability of cracking due 
to wind/current induced bending moments. The strength of the ice is such that it is very 
unlikely to yield to the relatively w r eak bending moments so induced. This conclusion is 
supported by the fact that significant noise levels have never been observed emanating 
from old inactive pressure ridge features. 



51 



4. Ice Fragment Bumping. 

This mechanism embraces all the bumping ancl grinding noises which would be 
expected to occur when ice fragments move about under an ice sheet. The probability 
that this mechanism is responsible for the creation of significant levels of ambient noise 
is reinforced by the recorded observation of audible sounds from the ice surface. Zubov 
(1943) quotes Weyprecht as stating that 

Sometimes shift at the bottom [ of the ridge keel ] is audible, with the complete re- 
pose of the ice at the top. This occurs probably as a result of the movement of water 
under the ice field. The difference of the movements of the ice field and of the water 
on which it lies, that is, the current of the water, is that power which levels the lower 
depth of the ice. 

Zubov (1943) also quotes Makarov on the existence and movement of ice blocks under 
the ice 

During the third Winter, a fissure was formed under the strip under the Fram. The 
fissure opened, lumps of significant size began to float out from below. "This 
shows", says Makarov, "that many lower lumps constantly travel. The current of the 
water and the movement of the ice change their direction so that the migrating ice 
block stops under certain conditions, under others it can move from the spot." 

It is neccessary to add that the lumps which compose the under-ice part of the 
hummocks [ ridges ] are not only transformed but are also gradually destroyed. 

Ice blocks or fragments are most likely to be the products of active ridging 
events. In the previous section the behaviour of young unconsolidated ridges was ex- 
amined quantitatively. The downward depression of material at the ridge axis was shown 
in Figures 21 and 22. This fractured ice appears to be thrust downward and trapped in 
buoyant equilibrium under the ice sheet by the ridging processes. Buck and Wilson 
(19S6) describe the forcing of ice blocks as large as four metres thick into the water 
column in their study of an active pressure ridge. 

It is reasonable then that given ice sheet motion relative to the water body be- 
low, the loose material will move horizontally. Indeed, further observations by Makarov, 
(Zubov, 1943) suggest that this is the process by which pressure ridge keels become 
vertically truncated and more widely spread in the horizontal than the corresponding sail 
formations. 



Given then that there are large ice fragments in motion along the rough under 
surface of an ice sheet, collisions, grinding and bumping are certain to occur. Noise from 
this source would probably be similar in nature to that observed as a result of floe-floe 
interactions in the marginal ice zone. Predominant frequencies are in the range 10 Hz to 
1 KHz according to Diachok and Winokur (1974). Certainly the impactive nature of the 
processes and the variety of fragment sizes is likely to yield noise spread over a wide 
frequency range. 

The intensity and perhaps the frequency range of ambient noise caused by this 
mechanism would be dependent on 

1. The relative current under the ice. 

2. The proliferation of ice debris under the ice sheet. 

3. The roughness of the under-ice surface. 

4. The thickness and density of the ice sheet. 

The current strength is of prime importance and is directly related to the translational 
speed of the ice and or oceanic or tidal currents. The presence of debris and ice 
roughness are both a function of the deformation history of the ice field and are likely 
to vary by region. Item 4 has been included as these factors are almost certain to have 
a bearing on the nature of the sound resulting from impactive events and its propagation 
into the sea. 

Modelling this mechanism is likely to be a complex task if approached from a 
purely physical standpoint. A more profitable approach might be to seek empirical re- 
lationships between the factors outlined above and measured ambient noise levels. This 
could be backed up by some simple field experiments involving perhaps acoustic moni- 
toring of ice blocks pushed through bore holes. 
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B. FLOW MECHANISMS. 

The potential noise creation mechanisms Falling into this category are all driven ex- 
clusively by the interaction of the sea with the underside of the ice sheet floating at the 
surface. The factors important in their operation are therefore the flow velocity of the 
water relative to the ice. the characteristics of the under-ice surface including cavity 
shapes and sizes, the roughness, and to some extent the thickness, of the ice cover. 

It seems intuiuively obvious that some noise, no matter how small, must be 
produced when water flows under rough ice. The object here is to explore the mech- 
anisms which may be responsible for this flow noise and assess their importance in the 
overall generation of ambient noise in an ice covered environment. 

Because of the complexity of the mechanisms outlined in this section, involving as 
they do, turbulent How over largely random surfaces, the prediction of the noise fields 
created is somewhat crude. Some of the physical processes are currently being researched 
and many of the concepts are still only poorly understood (Ross. 19S7). Where possible, 
estimates are made of likely sound characteristics in terms of frequency and intensity. 
These estimates are based on objective considerations but must be assumed accurate 
only to an order of magnitude. 

A typical composite Arctic ambient noise spectrum is presented in Figure 32. This 
serves to indicate the intensities required in order for a particular mechanism to be con- 
sidered important as a contributor. It is probable that How noise, whether or not it is 
found to be a significant component of the total noise under pack-ice, would be entirely 
masked by noise from other mechanisms at locations near the ice edge. For this reason, 
the measured spectrum chosen to represent the norm is from a central arctic pack-ice 
location. It is reproduced from Makris and Dyer (19S6) who made measurements at 
S3\Y20°£in April 19S2. 
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